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I. INTRODUCTION 

Lie's Group theory approach to the study of differen- 
tial equations has had a profound and enduring impact 
in Theoretical Physics, in particular in what respects to 
the invariance of differential equations under point trans- 
formations He showed that the one-dimensional free 
particle equation exhibits the symmetries corresponding 
to the eight-dimensional SL(3,R) group of point transfor- 
mations and that this is the maximum number of sym- 
metry generators for a second-order differential equation 
of the form [l[ 



V + h(y,y,x) = 0. 



(1) 



Of course, one often encounters situations in which one 
deals with a smaller number of generators. In a variety of 
such cases Q, equation {l) can nonetheless be recast in 
the form of a free particle equation by recourse to a non- 
local variable transformation. Now, as nonlocal transfor- 
mations can change both the number of symmetry gener- 
ators and the very physics of the original problem, they 
become a powerful tool for theoretic analysis. Einstein's 
field equations, whose nonlinear nature can indeed be 
faced with Lie's weaponry, constitute a particularly rele- 
vant example in this respect. 

It is via this equivalence (under nonlocal variable trans- 
formations) between different physical problem that we 
will search here for connections between 

• Cosmological problems in General Relativity, on 
the one hand, and 

• Information Theory Concepts pertaining to both 
Fisher's Information measure for translation fami- 
lies (FIM) [Hill and Tsallis' Nonextensive Ther- 
mostatistics (NET) [H-dl], on tne other one. 

This search is motivated by the recent findings of 
Ref. [25[ | : an intriguing relation between homogeneous 



and isotropic spatially flat cosmological models with 
no cosmological constant, on the one hand, and FIM- 
probability distributions that solve the diffusion equa- 
tion, as adapted a non- extensive (q = —1) Tsallis envi- 
ronment, on the other one. 

FIM's applications to diverse problems in theoretical 
physics have been pioneered by Frieden, Soffer, Nikolov, 
and others who have unveiled many FIM prop- 

erties and clarified its relation to Shannon's logarithmic 
information measure 



(2) 



where / is, of course, a probability density and < A > 
stands for the trace of (fA). 

General Relativity is a classical theory of gravitation 
that has been successfully applied for solving problems 
concerning astronomical scales of length and time, such 
as the formulation of cosmological models. Einstein's 
equation for the classical geometry adopts the appear- 
ance 



(3) 



where Rik is the Ricci tensor, R the scalar curvature and 
Tjft the energy-momentum tensor of matter and fields, 
with i,k = l,..,D, and D, the spacetime dimension. 

During the last few years much attention has been paid 
to models in which the universe, at very early stages of 
its evolution, was expanding 26-29]. The well known 
Big Bang model suggests a spatially homogeneous and 
isotropic space-time, stemming from a singularity in the 
remote past [3fjj . 

Here we will show that a generalization of the so-called 
Frieden-Soffer principle of Extreme Physical Information 



11], within the framework of a non extensive setting, 
may yield original insights into the workings of some of 
Einstein's field equations. 



2 



As a first step in such a direction, and in order to 
establish a convenient notation, we begin by refreshing, 
in Section II: i) the salient points of Fishers's measure 
tenets for translation families, ii) their extension to non- 
extensive settings, and iii) the Frieden Soffer Principle 
concerning extremalization of Physical Information [11] . 
In Section III we discuss the application of this Principle 
to a non-extensive environment. We present the main 
(cosmological) results of the present Communication in 
Section IV and, finally, some conclusions are drawn in 
Section V. 



II. A REVIEW OF FUNDAMENTAL NOTIONS 
A. The estimation theory bound 

In this work we shall concern ourselves exclusively with 
Fisher information measures I for translation families, 
that are invariant under Galilean and Lorentz transla- 
tions [3l[ . They refer to a measure of the inverse uncer- 
tainty in determining a position parameter by a maxi- 
mum likelihood estimation 32]. 

Estimation theory 5] provides one with a powerful re- 
sult that constitutes our stepping stone. Consider an 
isolated many-particle system that is specified by a po- 
sition parameter x and let f(x) describe the probability 
density (pd) for this parameter. The mean value of x 
for this pd is r\. Suppose the pd is unknown, and one 
wishes to determine it. Estimation theory asserts [j| that 
the best possible estimator rj est of our parameter, after 
a very large number of samples is examined, suffers a 
mean-square error e 2 from r\ that obeys a relationship 



involving Fisher's /, namely, 



V 



(4) 



where Fisher's information measure / (for tranlation fam- 
ilies) reads 



/ = 



(5) 



i.e. 



ill 

dx 



/O) 



(6) 



Any other estimator must have a larger mean-square er- 
ror. The only proviso to the above result is that all esti- 
mators r] est be unbiased, i.e., satisfy 



Thus, Fisher's information measure for translation 
families has a lower bound, in the sense that, no mat- 
ter what parameter of the system we choose to measure, 
/ has to be larger or equal than the inverse of the mean- 
square error associated with the concomitant measure- 
ment. This result, i.e., 



/ > 



1 



(8) 



is referred to as the Cramer-Rao bound, and constitutes 
a very powerful statistical result that was generalized in 
124] to a NET scenario. 



B. Estimation in a nonextensive Tsallis setting 

The phenomenal success of thermodynamics and sta- 
tistical physics crucially depends upon certain mathe- 
matical relationships involving energy and entropy, that 
can be translated in an essentially intact form into 
Tsallis' generalized Thermostatistics formalism [l5l - [l9j . 
whose main ingredient is the generalized entropy S q , that 
in terms of the real parameter q reads 



(1 - f- x ) ), qeTZ, 



(9) 



and coincides with Shannon's S for q — 1. Indeed, much 
work has been recently devoted to i) show that many of 
these relationships are valid for arbitrary values of Tsal- 
lis' parameter q, and ii) to find appropriate generaliza- 
tions for the rest. In this vein we just mention that, 
by suitably maximizing the entropy of the system, Cu- 
rado and Tsallis jl?} found that the whole mathemat- 
ical (Legendre-transform based) structure of thermody- 
namics becomes invariant under a change of the q— value 
(from unity to any other real number), while the con- 
nection of NET both with quantum mechanics and with 
Information Theory was established in [l8j , where it was 
shown that all of the conventional Jaynes-Boltzmann- 
Gibbs [33| results suitably generalize to the Tsallis' envi- 
ronment. For more details see [lSj . 

Of course, to verify that NET is useful, it is necessary 
to show that it appropriately describes certain physical 
systems with q- values that are different from unity. Much 
work in this respect has been performed recently (144 
refereed papers at the time of this writing) . We may cite 
applications to astrophysical problems flpL |20| , to Levy 
flights pllh to turbulence phenomena [22j . to simulated 
annealing [23j, etc. The interested reader is referred to 
for additional references. 

A suitable generalization of (0 to a NET setting 
should involve replacing any probability distribution / 
involved in evaluating expectation values by Tsallis' gen- 
eralized weights f q . As shown in [24j . one finds that 
Fisher's Information Measure in a Tsallis's nonextensive 



< Vest >= r). 



(7) setting (generalized FIM I q ) reads 
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(10) 



which can be abbreviate as GFIM. Instead of (J7J we have 
to deal, of course, with 

< Vest(y) > q = V- (11) 
In [|| one uses $TU\) to 

• generalize the Cramer- Rao bound @. 

• discuss connections between Tsallis' entropy, on the 
one hand, and Fisher's measure, on the other one. 

• derive a suitably generalized form of the Frieden- 
Nikolov's bound [HI to the time derivative 
This alternative form connects the entropy increase 

to Ig. 

C. The Frieden-Soffer EPI Principle 



The Principle of Extreme Physical Information (EPI) 
is an overall physical theory that is able to unify sev- 
eral sub-disciplines of Physics 11 1 . In Ref. [ll[ Frieden 



and Soffer (FS) show that the Lagrangians in Physics 
arise out of a mathematical game between an intelli- 
gent observer and Nature (that FS personalize in the 
appealing figure of a "demon" , reminiscent of the cele- 
brated Maxwell's one). The game's payoff introduces the 
EPI variational principle, which determines simultane- 
ously the Lagrangian and the physical ingredients of the 
concomitant scenario. 



FS 11| envision the following situation: some physical 
phenomenon is being investigated so as to gather suit- 
able, pertinent data. Measurements must be performed. 
Any measurement of physical parameters appropriate to 
the task at hand initiates a relay of information / (or I q ) 
from Nature (the demon) into the data. The observer 
acquires information, in this fashion, that is precisely I 
(or Ig). FS assume that this information is stored within 
the system (or inside the demon's mind). The demon's 
information is called, say, J (or J q ) [TTI | . 

Assume now that, due to the measuring process, the 
system is perturbed, which in turn induces a change 5 J 
(5J q )of the demon's mind. It is natural to ask ourselves 
how the data information / will be affected. Enters here 
FS's EPI: in its relay from the phenomenon to the data no 
loss of information should take place. The ensuing new 
Conservation Law states that SJ = SI, or, rephrasing it 



5(1 - J) = 0, 
so that, defining an action A 

A = I — J, 



(12) 



(13) 



EPI asserts that the whole process described above ex- 
tremalizes A. FS [ll[ conclude that the Lagrangian for 
a given physical environment is not just an ad-hoc con- 
struct that yields a suitable differential equation. It pos- 
sesses an intrinsic meaning. Its integral represents the 
physical information A for the physical scenario. On 
such a basis some of the most important equations of 
Physics can be derived [ll| . For an interesting Quantum 
Mechanical derivation see [34j . 

Within the present context our demon is working in a 
non-extensive fashion, so that the EPI principle should 
read 



An 



In 



Jo 



(14) 



i.e., using appropriate ^-generalizations of the intervening 
physical quantities. Thus the Frieden-Soffer information 
transfer game is played here according to 



5{I q - J q ) = 0. 



(15) 



III. THE WORKINGS OF FS'S DEMON IN A 
NON-EXTENSIVE SCENARIO 

A. Introductory remarks 

One starts here the FS game by using (fl"5"j). i.e., by 
extremalizing Fisher's generalized information for trans- 
lation families 



f( x y~ 2 f(x) 2 dx, 



(16) 



and considering a physical scenario in which the knowl- 
edge of the information demon [11] is limited to that 
concerning the normalization constraint. The demon's 



J q is 



Jo 



7o f{x) dx 7 



(17) 



where 70 is the pertinent Lagrange multiplier. The de- 
mon's information is reduced here to its minimum-core 
expression. 

Playing the Frieden-Soffer game then leads to 

- 2 f(x) + ( 2 -q)fPy- l0 f(x) 2 -" = 0, (18) 

g-dependent, non-linear differential equation that should 
yield our "optimal" probability distribution /. It easy to 
show that (|18p has the first integral 



/^ 2 ./ 2 + 7o/ = c 



and 



X — X{) = ± 



/ 



§-1 



:df 



(19) 



(20) 



V c - 70/ 

is the final expression for the general solution of (fT5|) . Of 
course, c and xq are arbitrary integration constants. 
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B. The q 7^ 1 environment 

^From (l20l) we see that c > 70/ in order to obtain a 
real probability distribution /. This condition holds in 
several circumstances, as follows: 

• For c > and 70 > 0, the probability distribution 
/ < c/70, i.e., it has an upper positive limit and it 
is bounded. 

• For c > and 70 < 0, the probability distribution 
/ > c/70, and thus it has a negative lower limit 
and it is unbounded. 

• For c < and 70 > 0, the probability distribution 
/ < c/70, i.e., it is negatively defined and we must 
reject it. 

• For c < and 70 < 0, the probability distribution 
/ > c/70, i.e., it has a positive lower limit and it is 
unbounded. 

We conclude that the first possibility is the only one 
that has a physical meaning. Indeed, inserting f(x m ) — 
c/70 into (|19|) . we obtain / = 0. Thus, at x rn the 
probability distribution / has an extremum. On the 
other hand, inserting these results into (fl~8|) we have 
/ = — -7t [if-] 9 (negatively defined). An /-maximum 
at x m ensues. 

The particular case c = is consistent with 70 < 
and needs separate analysis. Here, the momoparametric 
solution of (j2"0)l is given by 



Equation (|23j) arises in several interesting physical 
problems as, for instance, in the case of i) Einstein's 
field equations for homogeneous, isotropic and spatially 
flat cosmological models with no cosmological constant 
[29l [35M40l ] . and ii) Bianchi I- type metrics [4l|, for viscous 
fluids and a selfinteracting exponential potential scalar 
field. 

We shall tackle (|23l) by recourse to ideas presented by 
Chimento and Jakubi [351 ] . who show that its more gen- 
eral solution that can be explicitly writen as a function 
of x obtains for the special value a 2 //3 — 9. In our case, 
this assertion translates itself into two possible q- values, 
namely: 

• g = 2, and 

• q = -1, 

as will be discussed below. 

Following (35| , we perform the nonlocal change of vari- 
ables 



1 2 
- u 
2 



(26) 



dz 



and obtain 

drj 2 drj 

which is a linear ordinary differential equation, with 



+ (3 z 



0, 



(27) 



(21) 



drj — a u dx, 



(28) 



and 



This expression would yield a probability distribution 
that cannot be normalized if x ranges from minus to plus 
infinity. 



(29) 



The general solution of (f2"7| is of the form 



C. Relativistic scenario 



Evaluating derivatives (with respect x) in (fT5|i . and 
changing variables to 



z{ v ) = At e A+ " + A 2 e A " 



where 



(30) 



we obtain 



with 



and 



m 
m 1 



il + a u u + f3 u 3 = 



a = (2q-l) 



P = --q(l-q). 



(22) 

(23) 
(24) 

(25) 



2\a\ 



(31) 



are the roots of the pertinent characteristic equation, A\ , 
Ai being integration constants. 

Going all the way back from z to / we obtain the most 
general distribution function that verifies the Fricdcn- 
Soffer tenets. However, as far as one knows at the present 
stage (35l ], the road back can be traversed in explicit, 
analytical fashion only for 







a- 



9' 



(32) 
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which underlines the importance of the above mentioned 
special q- values (2 and —1) . In the next section we stress 
the connection between the probability distribution and 
the cosmological expansion factor, identifying both phys- 
ical quantities via equations (|23|) and ([27] ) . 



IV. COSMOLOGICAL APPLICATIONS 

Due to their nonlinear nature, exact solutions to Ein- 
stein's equations cannot easily be obtained. However, one 
finds diverse problems of physical interest where Einstein 
field equations (EFE) can be cast as particular instances 
of a second order nonlinear differential equation of the 
type 

V + ag(y)y + f3g(y) J dyg(y) + lg {y) = 0, (33) 

where y = y{x) and g{y) is a real function, a, /3, 7 are 
constant parameters. 

As examples of the preceding assertion, we may men- 
tion the case of EFE i) for homogeneous, isotropic and 
spatially fiat cosmological models with no cosmological 
constant [35l - [40| , ii) for a time decaying cosmological con- 
stant [42j . and iii) for Bianchi I- type metrics with a va- 
riety of matter sources 41 [ . 



A. The expansion rate of the Universe: different 
scenarios 

It is thought that quantum effects played a leading role 
in the early Universe, as, for instance, particle produc- 
tion and vacuum polarization, that arise within a quan- 
tal framework. It is known that these two phenomena 
can be modeled in terms of a classical bulk viscosity 
[43| . By recourse to the so-called Extended Irreversible 
Thermodynamics EH, a relativistic, second-order 
off-equilibrium approach, one can deal with a spatially 
flat homogeneous and isotropic space-time, filled with 
a causal viscous fluid and described by the Friedmann- 
Robertson- Walker line element 46] 



ds 1 



dt 2 - a 2 (t) [dx\ + dx\ + dx\\ , 



(34) 



with the full version of the transport equation for the 
bulk viscous pressure a 



T<J 



-3(H 




(35) 



In [35[ it was considered a causal viscous fluid whose 
equilibrium pressure obeys a 7-law equation of state, 
while the transport equation for the viscous pressure a 
was taken with e = 0, which corresponds to the so-called 
truncated theory. There, it was shown that H satisfies 
a modified Painleve-Ince equation that, in turn, can be 
cast into the form of (|3"3"|) . i.e., 



H + aHH + (3H 3 = 



(36) 



with f(y) = y and 7 = [2]. Comparing the latter with 
we get u = H, and, a fortiori a(t) = f{x) (from 
) . This result shows that the probability distribu- 
tion can be identified with the expansion parameter for 
the evolution of the universe by choosing the time vari- 
able t proportional to the information parameter x. This 
suggest that the solutions of (18) representing probabil- 
ity distributions without physical meaning, becuase they 
cannot be normalized, may be associated with the ex- 
panding (FRW) cosmological model. So that, we have a 
more extended background where the (FS)'s demon along 
with non-extensive scenario can be used. 

Cosmological models with a viscous fluid source have 
been also studied using the full causal Irreversible Ther- 
modynamics [H, |H, H3, S3- Relating the temperature 
T to the energy density in the simplest possible fashion 
that is able to guarantee a positive heat capacity, one 
finds that the expansion rate H verifies 



HH - (1 



AH H - BH = 



(37) 



where r > 0, while A, B and 70 are constants. By 
recourse to an adequate choice of variables (H = y n , 
with n = — 1/r) in Eq. (|37l) . the above equation becomes 
(again) identical to [2l.[29j. i.e., we can write 



ii + ay n y + Py 



2n+l 







(38) 



with f(y) — y n and 7 = 0. 

For a still different scenario consider now the early time 
evolution of a dissipative universe, in the relaxation dom- 
inated regime. One is again lead here to an equation for 
the expansion rate H that can be identified with (|2"3")l and 
(|3"3")) [48j |. Making once more the identification u = H, 
Eq.([3HD turns int0 E q- <[2"3"]>. so that a(t) = f(x), and 
the probability distribution can, as before, be identified 
with the expansion factor in the same manner as we did 
above. 

Let us discuss now the case of a perfect fluid source, 
satisfying an equation of state p = (a — l)p with cosmo- 
logical constant. The 00 Einstein equation 



where a{t) is the scale factor, H = a/ a is the expan- 
sion rate, r is the characteristic timescale for linear re- 
laxational effects, C is the linear bulk viscosity and T is 
the equilibrium temperature. 



3 a 



A 
I 



(39) 



transforms itself into Eq. (|33[) . with f(y) = y, and 
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B. Time variable and information variable 



H + (2 + 3a)HH + 3aH 3 - aAH = 



(40) 



provided we derive twice Eg. (1171) . Thus, a(t) — f(x) in 
a universe with no cosmological constant. 

Another interesting illustration of Eq. (|33|) can be ob- 
tained by considering an anisotropic Universe, described 
by a Bianchi type I metric, 



ds A 



J{t) 



-dt J + dz 1 



G(t) (« 



dx 1 



-pit) 



dy 



. ( 41 ) 

that is driven by a minimally coupled scalar field with an 
exponential potential Q. The Klein- Gordon equation for 
the scalar field and the Einstein equations for the metrics 
are expressed in terms of the semiconformal factor G and 
of their derivatives, as follows 



rn 
G 



kG 
2G' 



V= G< 



of - 



G 
2GV' 



G 1 G\ G , 1, 2 - 2 
G 2 \G } -G f+ 2 P = 



(42) 



(43) 



(44) 



(45) 



where m and a are arbitrary integration constants. The 
solutions for this set of equations can be obtained if one 
is able to solve the following equation for G 



C M 2 

G^ + {K-1)G+^- 
G G 



(46) 



where if is a constant parameter and M, c are integration 
constants. Once G(t) is known, one can compute the field 
cj>(t) and the remaining metric functions pit), f(t) from 
Eq. (|32|) . Eq. (|55J) . and Eq.([H]>. respectively. Making 



It should become clear by now that it is of great inter- 
est to investigate Eq. (|33[) from all conceivable points of 
view. Using nonlocal transformations this equation can 
be first linearized and then solved with the only proviso 
of the form invariance of Eq. (|33[) for an arbitrary function 

f{y). 

Now, if we relate y to the probability distribution func- 
tion / of the former Sections in the fashion 



y = 7 



and set 



H(t) = y, 



(49) 



(50) 



together with the choice: time variable t proportional to 
the information parameter x, one immediately realizes 
that equation ([3"3"l) translates itself into equation (|23p . the 
protagonist of our discussion in Section III. This shows 
that the probability distribution, in many circumstances, 
can be identified with the scale factor a(t) of the universe. 
In such a case, taking into account that the universe has 
a final, matter dominated Friedmann stage, the physical 
range of the Lagrange multiplier is given by 70 < 0. For 
c > the solutions stem from a singularity. However, 
for c < the solutions avoid the initial singularity. In 
the particular case c = the solutions (j2"Tj) represent a 
universe beginning or ending in a singularity at x = xq, 
for q > 1. For q < 1 the solutions diverge at the finite 
proper time x = xq- 

The special mapping: (information — > time) that 
arises in natural fashion from the present considerations 
allows one to give a degree of plasibility to the partic- 
ular Frieden-Soffer solutions with c = (f2"Tj) . since the 
scale factor, and its associated probability density, rep- 
resent the evolution of universe. In such circumstances, 
the scale factor is, of course, unbounded, it does not make 
any sense to try normalizing it. 

Summing up: we are leading to a variety of cosmolog- 
ical problems involving Einstein's field equations which 
via a nomlocal transformation and the special mapping: 
(information — > time) can be reinterpreted as a problem 
within a non-extensive scenario, by application of the 
nonextensive EPI Principle. 



G = y x ^ K , and r = —ct, 
Eq. (|4*o| becomes an equation of the type Eq. 



y" + y- n y> + ^y^ 
nci 



0. 



(47) 



(48) 



where a prime denotes the derivative with respect to r 
and n = l/K gj. 



V. CONCLUSIONS 

In the present communication we have achieved the 
following results 

• We translated the rules of the Frieden-Soffer Infor- 
mation Transfer game into a NET parlance. 

• By suitably playing the game, we found exact an- 
alytical solutions to special instances of Einstein's 
field equations. 
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• Nonlocal variable transformations of differential 
equations have allowed for a widening of the class 
of differential equations (particularly of the second 
order) which can be linearized and solved. 

In essence our procedure provides a method for deter- 
mining a useful nonlocal symmetry for a set of differ- 
ential equations characteristic of several different phys- 
ical problems. However, these problem can be seen to 
be are equivalent under nonlocal variable transforma- 
tions. In that peculiar sense, we see that some cosmo- 
logical solutions can be regarded as probability distri- 
butions, pertaining to a q =/= 1 NET environment, that 
maximize Fisher's information measure. One may think 
that Frieden-Soffer's demon is able to map the space-time 
continuum into a probability space. The solutions here 
studied can then be added to the impressive collection 
elaborated by Frieden and Softer [Til ]. 

That a NET environment may be the appropriate one 
in order to address physical problems connected with 
gravitation can be understood on the grounds that two 
systems A and B that interact in such a fashion can never 
be really "separated" (and become isolated ones) as ad- 
ditivity would demand, since they will always "feel" each 
other's presence on account of the infinite range of the 
associated force. Thus 

S(A+B) 7^ Sa + Sb (51) 

i.e., a non-extensive treatment becomes mandatory. In- 
deed, a classical problem related to Galactic models that 
had defied solution within an extensive Shannonian con- 
text has been successfully solved appealing to Tsallis' 
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NET [l9|. The present results should be viewed within 
such a background-landscape. 

As a last comment of this paper we mention that the 
results concerning to the equivalence of some expanding 
(FRW) cosmological model and the problem related with 
the (FS)'s demon in a non-extensive scenario are not con- 
sequence of the equivalence principle. They are a diret 
consequence of the nonlocal symmetries that have cer- 
tain nonlinear differential equations, in our case, the Ein- 
stein's field equations for a spatially flat (FRW) universe 
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